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Abstract 

Main theoretical uncertainties in estimating the indirect production of (6c)-quarkonium (B~ 
meson and its excited states) via top quark decays, t — > (be) + c + W + , are studied within the 
non-relativistic QCD framework. It is found that the dimensionless reduced decay width for a 
particular (6c)-quarkonium state, T n = T n /T t ^ w + +b , is very sensitive to the c-quark mass, while 
the uncertainties from the 6-quark and i-quark masses are small, where n stands for the eight 
(6c)-quarkonium states up to 0(v 4 ): {(bc^So)^, |(6c)( 3 ,Si)i), |(6c)( 1 Pi)i), |(6c)( 3 Pj)i) (with 
J = (1,2,3)), | (bc)^ So)sg) and \(bc)( 3 Si)8g) respectively. About 10 s tt-pairs shall be produced 
per year at CERN LHC, if adopting the assumption that all the higher Fock states decay to the 
ground state with 100% probability, then we shall have (l-038^Q 3 g2^ X 10 5 B~ events per year. 
So the indirect production provides another important way to study the properties of B~ meson 
in comparison to that of the direct hadronic production at CERN LHC. 
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I. INTRODUCTION 



The B c meson is a double heavy quark-antiquark bound state and carries flavors explicitly, 
which provides a good platform for a systematic studies of the 6 or c qua* decays. Since 
its first discovery at TEVATRON by CDF collaboration [1], B c physics is attracting more 
and wide interests. Many progresses have been made for the direct hadronic production 
of B c meson at high energy colliders g especiaHy, a computer program BCVEGPY for 
generating the B c events has been completed in Refs.[3|, m m an d has been accepted by 
several experimental groups to simulate the B c events. It has been estimated with the 
help of BCVEGPY that about 10 4 B c events are expected to be recorded during the first 
year of the CMS running with a lepton trigger [6j, and there are about 10 4 B c events with 
B c — > J/ty + 7T decays in three years of ATLAS running 

On the other hand, the indirect production of B+ or B~ , including its excited states, via t- 
decay or t-decay may also provide useful knowledge of these mesons. Without confusing and 
for simplifying the statements, later on we will not distinguish B^ and B~ (simply call them 
as B c ) and all results for B^t and B~ obtained in the paper are symmetric in the interchange 
from particle to anti-particle. With a predicted cross section for top quark pair production 
hundred times larger than at TEVATRON and a much higher designed luminosity, e.g. it is 
expected that at CERN LHC ~ 10 8 t t-pairs can be produced per year under the luminosity 
L = 10 34 cm~ 2 s~ x Q], the LHC is poised to become a "top factory". Therefore, the indirect 



production of B c throug 
properties of B c meson 



;h top quark decays shall provide another important way to study the 
9]. Within the non- relativistic QCD (NRQCD) framework [10], the 
decay channel t — > (be) + c+ W + in leading order (LO) a s calculation but with v 2 -expansion 
up to v 4 has been recently calculated with the so called 'new trace technology' [ll], where 
(6c)-quarkonium is in one of the eight Fock states: the six color-singlet states |(6c)( 1 S'o)i), 
{(bcj^Sili), K&c^Pi)].) and |(6c)( 3 Pj)!) (with J = (1,2,3)), and two color-octet states 
\(bc)( 1 So)sg) and | (&c) ( 3 S , i)s(?) respectively. It has been argued that when 10 8 ti events per 
year are produced at LHC, then it is possible to accumulate about 10 5 B c events per year 
via t-quark decays at LHC. Thus in comparison to that of the direct hadronic production, 
there may be some advantages in (6c)-quarkonium studies via the indirect production due 
to the fact that the top quark events shall always be recorded at LHC. 

Considering the forthcoming LHC running, and various experimental feasibility studies of 
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B c are in progress, it may be interesting to know the theoretical uncertainties quantitatively 
in estimating of B c production. The uncertainties of the direct hadronic production of B c 
through its dominant gluon-gluon fusion mechanism has been studied in Refs. [12I, rJ , while 
the present paper is served to study the uncertainties of the indirect mechanism through the 
decay channel t — > (be) + c + W + . These two cases are compensate to each other and may 
be useful for experimental studies. At the present, we shall restrict ourselves to examine 
the uncertainties at the lowest order, because the next-to-leading order (NLO) calculation 
cannot be available soon due to its complicatedness. For definiteness, we shall examine the 
main uncertainties that are caused by the value of the t-quark mass, the values of the bound 
state parameters m c and mj,, and the choice of the renormalization scale Q 2 . 

The remainder of the paper is organized as follows. Section II gives the calculation tech- 
nology for the indirect production of (6c)-quarkonium states through the top quark decays. 
Section III is devoted to present the numerical results and to discuss the corresponding un- 
certainties with the help of the formulae given in Ref. 11| . And section IV is reserved for a 
summary. 



II. CALCULATION TECHNOLOGY 

Within the non-relativistic QCD (NRQCD) frame work [10], the dimensionless reduced 
decay width for the production of (6c)-quarkonium through the channel t(po) ~~ > {bc)(pi) + 
C G°2) + W + (ps) takes the following factorization form: 



1 -if n (t-> (bc) + c + W + ) x { ° !!> 



(1) 



t^W++b N co i 

where T n stands for the reduced decay width for a particular (6c)-quarkonium state, and 
the sum is over all the (fec)-quarkonium states up to 0(v A ), which includes six color singlets 
KfcOOSo)!), l(^)( 3 Si)i>, |(fc)( 1 P 1 ) 1 ) and Kbc^Pj),) (with J = (1,2,3)), and two color 
octets \(bc)( 1 So)sg) and |(&c)( 3 Si) 8 g) respectively. N co i refers to the number of colors, n 
stands for the involved states of (6c)-quarkonium. N co i = 1 for singlets and N co i = N 2 — 1 
for octets. (O n ) stands for the decay matrix element that can be related with the wave 
function at zero -Rg(O) or the derivative of the radial wave function at origin R' P {0) through 
the saturation approximation [k]]. The overall factor l/T t ^ w + +b is introduced to cut off the 
uncertainty from the electroweak coupling. The decay width of the two body decay process 
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(bc)(pi) 



t(po) 




FIG. 1: Feynman diagrams for the indirect production of (6c)-quarkonium through top quark 
decays. 



t(Pi) ~ > b(p2) + W + (pz) that is dominant for the t-quark decays can be written as 

G F m 2 |p 2 | 



rv 



•W++b 



W2ir 



l-y 2 ) 2 + x 2 (l + y 2 -2x 2 



(2) 



where | p 2 1 = (1 — (x — y) 2 )(l — (x + y) 2 ), x = m w /m t and y = mb/m t . 

As shown in Fig.flT]), there are two Feynman diagrams for the concerned process t(po) — > 
(bc)(pi) + c(p 2 ) + W + (p 3 ). Due to the involved massive quarks, the calculation of the process 
is very complicated and lengthy, to simplify the calculation, we have improved a so called 



'new trace technology' to calculate the process [ll|]. Under such approach, we first arrange 
the whole amplitude into several orthogonal sub- amplitudes M ss i according to the spins of 
the t-quark (s') and c-quark (s), and then do the trace of the Dirac 7 matrix strings at 
the amplitude level by properly dealing with the massive spinors, which results in explicit 
series over some independent Lorentz-structures, and finally, we obtain the square of the 
amplitude. All the necessary formulae together with its subtle points for the square of the 
hard scattering amplitude H n (t — > (be) + c + W + ) can be found in Ref.jll], so we shall 



for more 



only present the main results here and the interesting reader may turn to Ref. 
detailed calculation technology. 

The involved color-singlet and color-octet matrix elements provide systematical errors 
for the NRQCD framework itself. Their values can be determined by global fitting of the 
experimental data or directly related to the wave functions at the zero point Rs{0) (or the 
derivative of the wave function at the zero point R' P (0)) derived from certain potential models 

nnnn 

for the color-singlet case, some potential models can be found in Ref. [14], [15|, [16|, UJ . A model 



dependent analysis of Rs(0) and R' P (0) can be found in Ref.[18|. where the spectrum of B c 

15J, the power-law potential 



under the Cornell potential [IJ], the Buchmiiller-Tye potential 



and the logarithmic potential [17| have been discussed respectively in their discussions, 
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ements 
we 



3, 



which shows that \R S (0)\ 2 G [1.508, 1.710]GeV 3 and \R' P {0)\ 2 G [0.201, 0.327]GeV 5 l . Since 
the model-dependent i?s(0) and -Rp(O) emerge as overall factors and their uncertainties can 
be conveniently discussed when we know their possible ranges well, so we shall not discuss 
such uncertainties in the present paper. More explicitly, we shall fix their values to be: 
\R S (0)\ 2 = 1.642 GeV 3 and [^(O)! 2 = 0.201 GeV 5 , which is derived under the Buchmiiller- 
Tye potential isj]. Secondly, although we do not know the exact values of the two decay 
color-octet matrix elements, (bc( 1 So)s\Os( 1 So) \bc( 1 So)s) and (bc^S^slOg^S^lbc^S^s), we 
know that they are one order in v 2 higher than the S'-wave color-singlet matrix e 
according to NRQCD scale rule. More specifically, based on the velocity scale rule 
have 

(bc( 1 S ) 8 \O 8 ( 1 S )\bcCS )s) ^ A s (v) 2 ■ (bcCSo^O^So^bcCSo)!) (3) 

and 

(bc( 3 S 1 ) 8 \Os( 3 S 1 )\bc( 3 S 1 ) 8 ) ~ A s (v) 2 • (6c( 8 5i) 1 |Oi( 8 5 1 )|6c( s 5 1 ) 1 > , (4) 

where the second equation comes from the vacuum-saturation approximation. As(v) is of 
order v 2 or so, and we take it to be within the region of 0.10 to 0.30, which is in consistent 
with the identification: A$(v) ~ a s (Mv) and has covered the possible variation due to the 
different ways to obtain the wave functions at the origin (S'-wave) and the first derivative of 
the wave functions at the origin (P-wave) etc. 

In addition to the color-singlet and color-octet matrix elements, the quark mass values 
m t , m c and also 'generate' uncertainties for the hadronic production. At present, these 
parameters cannot be completely fixed by fitting the available data of the heavy quarkonium. 
Furthermore, since the (fec)-quarkonium state is the non-relativistic and weak-binding bound 
state, we approximately have Mb c = vn^ + m c , which also is the requirement from the gauge 
invariance of the hard scattering amplitude. 

To choose the renormalization scale Q 2 is a tricky problem for the estimates of the LO 
pQCD calculation. If Q 2 is chosen properly, the results may be quite accurate. In the 
present case with three-body final state, there is ambiguity in choosing the renormalization 
scale Q 2 and various choices of Q 2 would generate quite different results. Such kind of am- 
biguity cannot be justified by the LO calculation itself, so we take it as the uncertainty of 



1 Since the Cornell potential has stronger singularity in spatially smaller states 18j . so we do not include 
its corresponding values for i?s(0) and R' P (Q). 
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TABLE I: Reduced decay width T n for the indirect production of B c through top quark decays 
with varying m c , rrt\, and wit, where n stands for a particular (c6)-quarkonium state. 



m c (GeV) 


1.5 


1.2 


1.8 


1.5 


1.5 


m b (GeV) 


4.9 


4.9 


4.5 


5.3 


4.9 


m t (GeV) 


172.5 


172.5 


172.5 


170. 


175. 


rosohtxio 4 ) 


3.590 


7.095 


2.053 


3.605 


3.573 


3.580 


3.600 


r ( a 5l)l (xio 4 ) 


4.975 


10.530 


2.690 


4.877 


5.065 


4.959 


4.991 


f (1Pl)l (xl0 5 ) 


3.543 


10.704 


1.449 


3.628 


3.470 


3.526 


3.559 


f (3Po)l (xl0 5 ) 


2.110 


5.294 


1.009 


2.301 


1.952 


2.106 


2.114 


f {3Pl)l (xl0 5 ) 


4.390 


12.528 


1.883 


4.585 


4.224 


4.372 


4.407 


f (3 p 2)l (xl0 5 ) 


4.718 


15.642 


1.765 


4.652 


4.774 


4.693 


4.743 



the LO calculation, although when the NLO calculation of the subprocess is available, the 
uncertainty will become under control a lot. While the NLO calculation is very complicated 
and it cannot be available in the foreseeable future, so here we take Q 2 as the possible char- 
acteristic momentum of the hard subprocess being squared. According to the factorization 
formulae, the running of a s should be of leading logarithm order, and the energy scale Q 2 
appearing in the calculation should be taken as one of the possible characteristic energy 
scales of the hard subprocess. As a default choice, we take Q 2 = Am 2 , since the intermediate 
gluon should be hard enough to produce a c and c pair as shown in Fig. ([I]). 

III. NUMERICAL RESULTS AND DISCUSSIONS 

Firstly, we study the uncertainties of m t , m c and mj, in £ a factorizable way' by fixing the 
renormalization scale Q 2 = Am 2 . For instance, when focussing on the uncertainties from m c , 
we let it be a basic 'input' parameter varying in a possible range m c = 1.5 ± 0.3 GeV with 
all the other factors, including the t-quark mass, 6-quark mass and etc. being fixed to their 
center values. The Particle Data Group value for the top quark mass is m t = 172.5±2.7 GeV 
[l9| . And the 6-quark mass m& varies within the region of m;, = 4.9 ±0.4 GeV. The reduced 
decay width f n for the indirect production of B c through top quark decays with varying m c , 
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m& and m t is shown in TABjTJ where n stands for a particular color singlet (cfe)-quarkonium 
state. The results for the two S'-wave color octet can be conveniently obtained from that 
of color singlet S-wave (c6)-quarkonium states and by setting A s {y) G [0.10,0.30]. The 
second column of TABjT]is for the center values of all these parameters, the third and fourth 
columns setting the upper and the lower limit for m c varying within the region of [1.2, 1.8] 
GeV, the fifth and the sixth columns setting the upper and the lower limit for varying 
within the region of [4.5, 5.3] GeV, and the seventh and eighth columns setting the upper 
and the lower limit for m t varying within the region of [170., 175.] GeV respectively. 

From TABjlj it is found that the reduced decay width f n is very sensitive to m c . f n 
decreases with the increment of m c , and more definitely, when m c increase by steps of 
O.lGeV, F n decreases by 10% ~ 20% for S'-wave states and by 25% ~ 35% for P-wave states. 



This condition is similar to the direct hadronic production [12J, which is caused by the fact 
that a larger m c leads to a smaller allowed phase space. Summing up all the mentioned 
Fock states' contribution, we obtain Enf„ = (l.038tjj;|») x 10 ~ 3 for m c e [1.2, 1.8] GeV 
and Ag(i>) G [0.10,0.30], where the center value is for m c = 1.5GeV, m b = 4.9GeV, m t = 
172.5GeV and A s (v) = 0.2. 

The reduced decay width f n slightly decreases with the increment of mj for n = 
K&cX 1 ^)!), n = \(bc)( 1 S ) 8 g), \(bc){ 1 P l ) 1 ), |( 3 P )i) and {{bc^P^) respectively, but in- 
creases with the increment of mb for n = |(6c)( 3 S'i)i), | (6c) ( 3 S , i)s(?) and |(frc)( 3 P2)i) re- 
spectively. As for the direct hadronic production of B c , since there is a 6-quark jet in the 



final state, so the production shall always decrease with the increment of 6-quark mass [12]. 
While for the present case, there is no such 6-quark jet in the final state, so the condition is 
slightly different. Further more, it is found that when increase by steps of 0.2GeV, the 
uncertainties is less than 1%. Summing up all the mentioned Fock states' contribution, we 
obtain En f „ = (l.038±g;^) x 10~ 3 for m b G [4.5,5.3]GeV and A s (v) G [0.10,0.30]. 

The reduced decay width f n slightly increases with the increment of m t , which is due to 
the , arger phase spaee for a large, m t . To check the results of Ref.Q, we alsocalculate the 
results for m t = 176GeV, which shows a good agreement with those of Ref . [1 ll] . Summing 
up all the mentioned Fock states' contribution, we obtain Z) n f n = (l-038lo;o2g) x 10 -3 for 
mt G [170, 175]GeV and A s {v) G [0.1,0.3]. 

Secondly, we study the uncertainties caused by the various choices of Q 2 , where for consis- 
tency, the leading order a s running is adopted, i.e. a s {Q 2 ) = 47r/[(ll — |n/) \n(Q 2 / Aq CD )], 
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TABLE II: Reduced decay width f n for the indirect production of B c through top quark decays 
with three typical renormalization scale Q 2 , where n stands for a particular color-singlet (cb)- 
quarkonium state. 



m c (GeV) 


Q 2 = E 2 Bc 


Q 2 = 4mg 


Q 2 = Am 2 c 


r (1 s o)l (xio 4 ) 


2.125 


1.738 


3.590 


f (35l)l (xl0 4 ) 


2.875 


2.408 


4.975 


r (1Pl)l (xl0 5 ) 


2.040 


1.715 


3.543 


f (3Po)l (xl0 5 ) 


1.177 


1.022 


2.110 


f (3Pl)l (xl0 5 ) 


2.546 


2.125 


4.390 


f ( 3 P2)l (xl0 5 ) 


2.700 


2.284 


4.718 




FIG. 2: Differential distributions dT /ds\ and dT/ds2- The shaded shows the uncertainty and the 
solid line is for the center value with m c = 1.5GeV, nib = 4.9GeV and rrit = 172.5GeV. 

where rif — 3 and Aqcd = 200MeV. We choose three typical renormalization scale Q 2 : 
Type A: Q 2 = 4m 2 c ; Type B: Q 2 = 4m 2 b ; Type C: Q 2 = E% c , where E Bc stands for the B c 
meson energy in the top quark rest frame, and by setting S2 = (p2 + P3) 2 for the channel 
t(p ) — > (bc)(pi) +c(p 2 ) + W + (p 3 ), we have E Bc = (m 2 + M Bc — s 2 )/(2m t ). The uncertainties 
for the reduced decay width f n with three typical choices of Q 2 are given in TAB Jill where 
m c = 1.5GeV, m b = 4.9GeV and m t = 172.5GeV. It is found that the reduced width for 
Q 2 = Ami i s only about half of that of Q 2 = ^m 2 c) which is a comparatively large effect. 
Finally, we discuss the combined effects of all the above mentioned uncertainty sources 
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FIG. 3: Differential distributions dT / dcos6\2 and dT /dcosO^. The shaded shows the uncertainty 
and the solid line is for the center value with m c = 1.5GeV, m& = 4.9GeV and mt = 172.5GeV. 



by varying m c G [1.2, 1.8] GeV, mj G [4.5, 5.3]GeV, m t G [170, 175]GeV and by taking one 
of the three typical renormalization scales simultaneously, where all the mentioned Fock 
states' contributions shall be summed up. Additionally, the value of m& + m c can not be 
too small, as has been found both experimentall y a nd theoretically that the mass of the 



ground state (fec)-quarkonium is around 6.30GeV [20j, |2l(, so we imply m& + m c ^ 6.20GeV 



as an extra constraints. Summing up all the mentioned Fock states' contribution, we obtain 
J2 n ^n = (l-038lo!78l) x 10~ 3 . Let us show some more characteristics of the decay t — > 
(be) + c + W + . The differential distributions of the reduced decay width versus the invariant 
masses si = (pi + P2) 2 and S2 = (P2 + P3) 2 , i-e. dT/dsi and dT/ds2 are shown in Fig. ([2D - 
While the differential distributions of the reduced decay width versus cos #13 and cos#i 2 , i.e. 
dT /d cos 9 12 and dV/d cos 6*13 are shown in Fig. where 6*13 is the angle between pi and 
p s , and 9 12 is the angle between pi and p 2 respectively in the t-quark rest frame (p = 0). 
In drawing the curves, all the mentioned Fock states' contribution have been summed up 
for convenience. The shaded band shows the corresponding uncertainty, with the upper 
edge of the band is obtained by setting m c = 1.2GeV, nib = 5.0GeV, m t = 175GeV and 
Q 2 = Ami and the lower edge of the band is obtained by setting m c = 1.8GeV, m& = 5.3GeV, 
m t = 170GeV and Q 2 = Am 2 , and the center solid line is for m c = 1.5GeV, m& = 4.9GeV, 
m t = 172.5GeV and Q 2 = Am 2 c . 
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IV. SUMMARY 



In the paper we have presented quantitative studies on the main uncertainties in estimat- 
ing the indirect production of the (6c)-quarkonium via top quark decays, t — > (be) + c + W + . 
It is found that the reduced decay width f n is very sensitive to the c-quark mass, while 
the uncertainty from the 6-quark and t-quark masses are small. The renormalization scale 
also affects the decay width to a certain degree. A comparative study on the similarity and 
difference of the direct and indirect production has also been presented in due places. About 
10 8 tt-pairs shall be produced per year at CERN LHC, if adopting the assumption that all 
the higher Fock states decay to the ground state with 100% probability, then we may have 
(l.038^;fl) x 10 5 B-{B+) events per year. So the indirect production is another impor- 
tant way to study the properties of B c meson in comparison to that of the direct hadronic 
production at LHC. Further more, the contribution from the P-wave states together with 
the two color-octet Fock states can be about 20% in total, so the P-wave production itself 
is worthwhile to study the possibility of directly measuring the P-wave B c states. 
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